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Abstract 

A new iterative algorithm for solving initial data inverse problems from par- 
tial observations has been recently proposed in Ramdani, Tucsnak and Weiss |15j . 
Based on the concept of observers (also called Luenberger observers), this algo- 
rithm covers a large class of abstract evolution PDE's. In this paper, we are con- 
t-h cerned with the convergence analysis of this algorithm. More precisely, we provide 

a complete numerical analysis for semi-discrete (in space) and fully discrete approx- 
imations derived using finite elements in space and finite differences in time. The 
analysis is carried out for abstract Schrodinger and wave conservative systems with 

^ bounded observation (locally distributed). 
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1 Introduction 

The goal of this paper is to present a convergence analysis for the iterative algorithm 
recently proposed in Ramdani, Tucsnak and Weiss [Tb\ for solving initial state inverse 
problems from measurements over a time interval. This algorithm is based on the use 
back and forth in time of observers (sometimes called Luenberger observers or Kalman 
observers; see for instance Curtain and Zwart [3]). Let us emphasize that during the 
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last decade, observers have been designed for linear and nonlinear infinite-dimensional 
systems in many works, among which we can mention for instance Auroux and Blum [1] 
in the context of data assimilation, Deguenon, Sallet and Xu [5], Guo and Guo [H], Guo 
and Shao [9] in the context of wave-type systems, Lasiecka and Triggiani [1 2j , Smyshlyaev 
and Krstic [17] for parabolic systems and Krstic, Magnis and Vazquez [10J for the non 
linear viscous Burgers equation. 

Let us first briefly describe the principle of the reconstruction method proposed in [T5] 
in the simplified context of skew-adjoint generators and bounded observation operator. 
We will always work under these assumptions throughout the paper. Given two Hilbert 
spaces X and Y (called state and output spaces respectively), let A : V (A) — > X be 
skew-adjoint operator generating a Co-group T of isometries on X and let C e C(X, Y) 
be a bounded observation operator. Consider the infinite dimensional linear system given 
by 

/ z(t) = Az(t), Vt>0, 

\y(t) = Cz(t), VtG[0,r]. ^ 

where z is the state and y the output function (throughout the paper, the dot symbol is 
used to denote the time derivative). Such systems are often used as models of vibrating 
systems (e.g., the wave equation, the beam equation,...), electromagnetic phenomena 
(Maxwell's equations) or in quantum mechanics (Schrodinger's equation). 




t = T 



t = 



Figure 1: An initial data inverse problem for evolution PDE's : How to reconstruct the 
initial state (light grey) for a PDE set on a domain Q from partial observation on O x [0, r] 
(dark grey)? 

The inverse problem considered here is to reconstruct the initial state z = z(0) of 

Fig. 0. 



system (1.1) knowing (the observation) y(t) on the time interval [0,r] (see Fig. yy). Such 
inverse problems arise in many applications, like thermoacoustic tomography Kuchment 
and Kunyansky [H] or data assimilation Puel [H]. To solve this inverse problem, we 
assume here that it is well-posed, i.e. that (A, C) is exactly observable in time r > 0, i.e. 
that there exists k T > such that 



\y(t)\\ 2 dt> k 2 T \\z 



Vz e V(A). 
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Following Liu [TBI Theorem 2.3.], we know that A + = A — C*C (respectively A~ = 
—A — C*C) generate an exponentially stable Co-semigroup T + (respectively T~) on X. 
Then, we introduce the following initial and final Cauchy problems, called respectively 
forward and backward observers of (1.1) 

z+{t) = A+z + (t) + C*y(t), Vt G [0, r], . . 

z+(0) = 0, 1 } 

z-(t) = -A-z-(t)-C*y(t), Wg[0,t], 
z~(t) = z + (t). 

Note that the states z + and z~ of the forward and backward observers are completely 
determined by the knowledge of the output y. If we set L r = T~T+, then by [TS1 
Proposition 3.7], we have rj := ||L T ||,c(x) < 1 and by [15, Proposition 3.3], the following 
remarkable relation holds true 

z = (I-h T )- 1 z-(0). (1.4) 

In particular, one can invert the operator (J — L T ) using a Neumann series and get the 
following expression for the initial state 

2o = X)L^-(0). (1.5) 

n=0 

Thus, at least theoretically, the reconstruction of the initial state is given by the above 
formula. Note that the computation of each term in the above sum requires to solve the 
two non-homogeneous systems (1.2) and (1.3). In practice, the reconstruction procedure 
requires the discretization of these two systems and the truncation of the infinite sum 
in (1.5) to keep only a finite number of back and forth iterations. For instance, if we 
consider a space semi-discretization corresponding to a mesh size h (typically a finite 
element approximation), one can only compute 

N h 

Z ,h = J2 L lr^(0), (1.6) 

n=0 

where 

• hh, T = T^ T T^ T , where T^ r G C(X) are suitable space discretizations of T^, 

• z^ (0) G X h is an approximation of z~(0) in a suitable finite dimensional subspace 
X h of X, 

• Nh is a suitable truncation parameter. 

Similarly, if a full discretization described by a mesh size h and a time step At is consid- 
ered, one can compute 

Nh.At 
n=0 

where 
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At,K — At K^t At ki where T^ At K are suitable space and time discretizations 



of T d T , 



T 5 
,0 



• (zfr) G X h is an approximation of z~(0), 

• Nh,At is a suitable truncation parameter. 

For the sake of clarity, the precise definitions of the spaces and discretizations used will 
be given later in the paper. 

Our objective in this work is to propose a convergence analysis of z 0t h and Zo,/i,At 
towards Zq. A particular attention will be devoted to the optimal choice of the truncation 
parameters and iV^ a* f° r given discretization parameters (mesh size h and time step 



At). Let us emphasize that our error estimates (see (2.8), (2.27), (3.15) and (3.36)) 
provide in particular an upper bound for the maximum admissible noise under which 
convergence of the algorithm is guaranteed. As usually in approximation error theory of 
PDE's, some regularity assumptions are needed to obtain our error estimates. Namely, 
our result allows us to reconstruct only initial data contained in some subspace of X 
(namely T>(A 2 )). Moreover, our analysis only holds for locally distributed observation 
(leading to bounded observation operators). 

Throughout the paper, we denote by M a constant independent of r, of the initial 
state zq and of the discretization parameters h and At, but which may differ from line to 
line in the computations. 

The paper is organized as follows: in Section [2] we provide a convergence analysis 
of the algorithm for an abstract Schrodinger type system, by considering successively 



the semi-discretization (Subsection 2.1) and the full discretization (Subsection |2.2). In 



Section |3j similar results are given for an abstract wave system. Once again, we tackle 



successively the semi-discretization (Subsection 3.1) and the full discretization (Subsec- 
tion 3.2). Finally, the Appendix is devoted to the proof of two technical lemmas which 



are used several times troughout the paper. 



2 Schrodinger equation 

Let X be a Hilbert space endowed with the inner product (•, •). Let A : T> (A ) — y X be 
a strictly positive self-adjoint operator and C G C(X, Y) a bounded observation operator, 
where Y is an other Hilbert space. The norm in T>(Aq) will be denoted by || • || a . We 
assume that there exists some r > such that (iAq, C) is exactly observable in time 
r. Thus by Liu [131 Theorem 2.3.], A+ = iA - C*C (resp. A~ = -iA - C*C) is the 
generator of an exponentially stable Co-semigroup T + (resp. T - ). We want to reconstruct 
the initial value z of the following system 

/ z{t) = iA z{t), Vt^O, 

\y(t) = Cz(t), VtG[0,r]. {lA} 

Throughout this section we always assume that zq G T> (Aq). Thus by applying Theorem 
4.1.6 of Tucsnak and Weiss [18], we have 

zeC([0,r],V(Al))nC l ([0,r],V(A )). 
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The forward and backward observers (1.2) and (1.3) read then as follows 



z+(t) = iA z+(t) - C*Cz + {t) + C*y{t), Vt G [0, r] 
z+(0)=0, 



(2.2) 
(2.3) 



z-(t) = iA z-(t) + C*Cz~(t) - C*y(t), Vt G [0,r], 
z~(t) = z + {t). 

Clearly, the above systems can be rewritten in the general form of an initial value Cauchy 
problem (simply by using a time reversal for the second system) 



g(t) = ±iA q(t) - C*Cq{t) + F{t), Vt G [0, r], 
9(0) = 9o, 



(2.4) 



where we have set 



for the forward observer g2j) : F(t) = C*y(t) = C*Cz(t) and q = 0, 

for the backward observer ( |2.3[ ) : F(t) = C*y(r — t) = C*C,2(t - £) and q 
z+(r)eV(Al). 



2.1 Space Semi-Discretization 
2.1.1 Statement of the main result 



We use a Galerkin method to approximate system (2.4). More precisely, consider a family 
(X h ) h>0 of finite-dimensional subspaces of V (Aq J endowed with the norm in X. We 



denote tt/j the orthogonal projection from X> |^4q J onto X^. We assume that there exist 
M > 0, 9 > and h* > such that we have for all ft, G (0, ft,*) 



\ith<p — vll ^ Mft. 



W>GD L4 2 . 



(2.5) 



Given go G P(Ag), the variational formulation of (2.4) reads for all t G [0,r] and all 
(/? G P A 2 as follows 



(g(t), ^) = ±i (q(t), <p)± - (C*Cq(t), <p) + (F(t), <p) , 
9(0) = go- 



(2.6) 



Suppose that go,h G -A^ and Fh are given approximations of go and F respectively in the 
spaces X and L 1 ([0, r], X). For all t G [0, r], we define qh(t) G X^ as the unique solution 
of the variational problem 



(qh(t),<Ph) = ±i (qh{t),Vh)± ~ (C*Cq h {t),ip h ) + (F h {t),ip h ) 
Qh(0) = Qo,h- 



(2.7) 



for all ip h e X h . 

The above approximation procedure leads in particular to the definition of the semi- 
discretized versions T± of the semigroups that we will use. Indeed, we simply set 



T f + 9o ^ T+ t g = q h {t) 



T t 9o ^ T h t q = q h (r- t) 



where is the solution of equation (2.7) with the corresponding sign and for F h = and 
Qo,h — 7T/i?o- The approximation of L T = T^T^ follows immediately by setting 

W = TJ^ T T+ T . 

Assume that is an approximation of the output y in -^ 1 ([0, r], F) and let ^ and 
denote the Galerkin approximations of the solutions of systems (2.2) and (2.3), satisfying 
for all t G [0, r] and all G Xh 

(zUt),<p h ) = i (z+(t), <p h ) h - (C*Cz + (t), Vh ) + (C*y h (t),<p h ) , 
4(0) = 0. 

(z^(t),(p h ) =i{zZ(t),<p h )i + (C*Cz h (t), Vh ) - (C*y h (t),<p h ) , 



Thus, our main result in this subsection reads as follows. 

Theorem 2.1. Let Aq : T> (Aq) — )■ X be a strictly positive self-adjoint operator and 
C G C(X,Y) such that C*C G C(V(A 2 )) nC(V(A )). Assume that the pair (iA ,C) 
is exactly observable in time r > and set rj := ||L r ||£nn < 1. Let zq G T> (Aq) be the 
initial value of (2.1) and z 0) h be defined by (1.6). 

Then there exist M > and h* > such that for all h G (0, h*) 



\z - z h\\ < M 



n N h +l 

Y^ + h e rN 2 h ) j|:„|| 2 ■■ ■ X h 



\C* (y(s)-y h (s))\\ds 



A particular choice of Nh leads to an explicit error estimate (with respect to h) as 
shown in the next Corollary (the proof is left to the reader because of its simplicity) 



Corollary 2.2. Under the assumptions of Theorem \2. 1\ we set 



In rj 



Then, there exist M T > and h* > such that for all h G (0, h*) 



Zo-Zo, h \\<M T [h d \n 2 h\\z \\ 2 + \lnh\ / \\C* (y(s) - y h (s)) \\ds 



(2.8) 



Remark 2.1. In fact, Theorem 



2.1 



still holds true for z G V I^Aq J (with the same proofs 
Nevertheless, we have not been able to carry out this 



and slightly adapting the spaces) - 
analysis for the fully discrete approximation in this case. This is why we restricted our 
analysis to the case of an initial data zo G T> (Aq). 
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2.1.2 Proof of Theorem 1231 



Before proving Theorem |2.1[ we first need to prove some auxiliary results. The next 
Proposition, which constitutes one of the main ingredients of the proof, provides the 



error estimate for the approximation in space of the initial value problem (2.6) by using 



the Galerkin scheme (2.7). 



Proposition 2.3. Given qo G T> (Aq) and go,/i E Xh, let q and qh be the solutions of (2.6) 



and (2.7) respectively. Assume that C*C G C(V(A )). Then, there exist M > and 



h* > such that for all t G [0, r] and all h G (0, h*) 



h h q{t) ~ q h {t)\\ < \\n h q - q Qh \\ + Mh e 



*(ll?o| 



|-^||l,oo) + £ 2 ||-F||2,oo 



\F{s)-F h {s)\\ds. 



Proof. First, we substract (2.7) from (2.6) and obtain (we omit the time dependence for 
the sake of clarity) for all tph G Xh 

(q - q h ,(Ph) = ±* (<? ~ qh,<Ph)i - (C*C(q - q h ),Vh) + (F - F h ,ip h ) . 

Noting that (nhq — q, <ph) i = for all (fh G Xh and that •n^q makes sense by the regularity 



of q (see (4.1)), we obtain from the above equality that for all (fh G Xh 



(n h q ~ qh, <Ph) = {^hq ~ q, V>h) ± « (^hQ - qh, fh) i 

-(CTC(q-q h ),<p h ) + {F-F h ,<p h ). (2.9) 

On the other hand, setting 



we have 



£h = \ hhq ~ qhf, 



t h = Re (ir h q - q h , rc h q - qh) ■ 



Applying (2.9) with iph = TChq — qh and substituting the result in the above relation, we 



obtain by using Cauchy-Schwarz inequality and the boundedness of C that there exists 
M > such that 

£h < (\\TT h q ~ q\\ + M\\ir h q - q\\ + \\F - F h \\) \\ir h q - qh\[ - 
Sf d 

Since = —y/2Sh, the integration of the above inequality from to t yields 

y/2£ h dt 



K<?(0 ~ qh (t)\\ < W^hQo - Qo,h\\ + / ihhq(s) - q(s)\\ + M\\7r h q(s) - q(s)\\) ds 

Jo 

+ / \\F(s) - F h (s)\\ds. (2.10) 
Jo 
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Thus, it remains to bound \\iT h q(t) — q(t)\\ and ||7r^g(t) — q(t)\\ for all t G [0, r]. Using 



(2.5) and the classical continuous embedding from V(A a ) to T)(A") for a > (3, we get 



that 



\n h q(t) - q(t)\\ < Mh 6 \\q(t)\\ i < Mh $ \\q(t)\\!, 
\7T h q(t) - q(t)\\ < Mh°\\q(t)\\l < Mh%(t)\\ 2 , 



Vt G [0,r], /i G (0,/i*). 



Using relations (4.2) and (4.3) proved in Lemma 4.2 of the Appendix, we get for all 



t G [0,r] and all h G (0,/i*) 

IKg(t) -g(t)|| + ||7r h g(*)-g(0|| < Mh° (\\q \\ 2 + t\\F\\ 2j00 + \\F\\ 1;00 ) . 
Substituting the above inequality in ( |2.10[ ), we get the result. ■ 

Using the last result, we derive an error approximation for the semigroups T 1 * 1 and for 
the operator L t = T^Tf. 



Proposition 2.4. Under the assumptions of Proposition 2.3 , the following assertions 
hold true 

1. There exist M > and h* > such that for all t G (0, r) and all h G (0, h*) 

\\7T h T+q - T+ t q \\ < Mth e \\q \\ 2 . (2.11) 

||7r At T-g -T M g || < M(r — t)h 9 \\q \\ 2 . (2.12) 

2. There exist M > and h* > such that for all n G N, all t G [0, r] and all 
h G (0, h*), we have 



\\L r ; qo - L£ tg o|| < M(l + nr)/i || go ||2. 

Proof. 

1. It suffices to take F = F ft = and go.ft = t^Mo i n Proposition 2.3 

2. We first note that 



|L^ - KM\ ^ W h tQo ~ TThKQoW + \WhKQa ~ K,tQol 



(2.13) 



(2.14) 



Using (2.5) and the fact that ||lLt||£(x>(A)) — 1 proved in Lemma 4.1 of the Appendix, the 
first term in the above relation can be estimated as follows 



\L?q - n^go\\ < Mh e \\q \\ 2 , V/i G (0, h*). 



For the second term in (2.14), we prove by induction that for all n G N 



IKL^d ~KM\ ^ Mnrh d \\q \\ 2 , Vh G (0,/i*). 

By definition, we have 

\\ir h h t q - h hit q \\ = \\n h TtTtq - T^ t T+ t g ||, 

< ||7r h T t -T+g - T^T+g || + ||T^(T+g - T+ t9o ) 



(2.15) 



(2.16) 



S 



By Lemma 4.1 of the Appendix and equation (2.12), we get 



KT-T+go - T^T+goll <M(r- t)h e \\q \\ 2 , V/i G (0, h*). 

Obviously ||Tj~||£nn is uniformly bounded with respect to h (this follows for example 
from (2.12)), and thus by (2.5) and equation (2.11), we have 

||T M (T+g - T+ g )|| < ||T+g - 7T fe T+g || + KT+g - T+ g || 

< Mth e \\q \\ 2 , \/he{0,h*). 



Consequently 



\7ihLtqo - h h>t q \\ < MTh e \\q \\ 2 , Wh G (0, h*), 



(2.17) 



(2.18) 



which shows that (2.16) holds for n — 1. Suppose now that for a given n > 2, there holds 

KLrto - L^-^oll < M{n - l)rh 6 \\q \\ 2 . 

We write 

||7r^g - g || < H^W^o - L^L^goll + WU^t'^o - K^Vo) 



Thanks to Lemma 4.1 and to the uniform boundedness of HL^H^x) with respect to h 
(which follows from the uniform boundedness of ||T^J|) and using (2.17) and (2.18), we 
obtain 



\n h L,?q - KM < M(t + (n- l)r)h 9 \\q \ 



2- 



which is exactly (2.16). Substituting (2.15) and (2.16) in (2.14), we obtain the result. 



We are now able to prove Theorem 2.1 



of Theorem 2J_, Introducing the term K,t z (0)> we rewrite z — z ,/i in the following 
form 



n=0 



oo N h 

z -z 0)h =^L>-(0)-^L^(0), 

n=0 n=0 

N h N h 

= L >"(o) + E ( L " - L ir) z-(o) + E L ^ (*"(°) - ** (°)) • 

n>N h 



n=0 



n=0 



Therefore, we have 
where we have set 



\zo - Z h\\ < Si + 5*2 + 5*3 



(2.19) 



' E ll L >~(°)ll> 

n>N h 
N h 

5 2 = EIK L "- L U ^ co) 1 1 



n=0 

N h 



EII L U 



v n=0 



C(X) 

9 



|*-(0)-*fc(0)| 



Note that the term S\ is the truncation error of the tail of the infinite sum (1.5), the term 
S 2 represents the cumulated error due to the approximation of the semigroups T while 
the term 5*3 comes from the approximation of the first iterate z~(0) of the algorithm. 

Since r\ = ||L r ||£(x) < 1, using relation (1.4), the first term can be estimated very 
easily 



Si < M 



V 



N h +1 



7] 



\ Z 2- 



(2.20) 



The term S 2 can be estimated using the estimate (2.13) from Proposition 2.4 



S 2 < M [5^(1 + nr) J ^||z-(0)|| 2 , \/h G (0, h*). 

\n=0 / 



Therefore, using (1.4) and the fact that ||L r ||x)(A 2 ) < 1 (see Lemma 4.1) in the above 
relation, we finally get that 



S 9 . < M 



1 + (1 + r)N h + N 2 h r h e \\z \\ 2 , V/i G (0, h*). 



(2.21) 



It remains to estimate the term S 3 . As 77 = ||L r ||£(x) < lj (2.13) implies that ||L/j T ||£(x) 
is also uniformly with respect to h bounded by 1, provided h is small enough. Hence, we 

U<U< ■% < -UA- /( ||--(0)- ^(0)1, (2 22) 



< MN h (\\z-(0) - n h z-(0)\\ + \\7r h z-(0) - ^(0)||) . 



By using (2.5) and (1.4), we immediately obtain that 

|U~(0) - 7r^ - (0)|| < Mh e \\z \\ 2 . 



(2.23) 



To estimate the second term -n^z (0) — z h (0), we apply twice Proposition 2.3 first for 
the time reversed backward observer z~{r — •) and then for the forward observer z + (the 



time reversal step is introduced as in the formulation of Proposition |2.3[ only initial value 
Cauchy problems can be considered). After straightforward calculation we obtain that 
for all h G (0,/i*) 



k & 2-(0)-^(0)|| <Mh 6 r(||^ + (r)|| 2 +||C*y|| li0O )+r 2 ||C*|/| 



2,00 



+ 



\C*{y{r-s)-y h (T-s))\\ds+ \\C* (y(s) - y h (s)) \\ds. (2.24) 



Applying (4.2) of Lemma 4.2 of the Appendix with zero initial data, we obtain that 

lk + (r)||2<r||C*y|| 2iCX) . 



Therefore (2.24) also reads 

\\n h z-(0)-z h (0)\\<Mh\r + T 2 )\\C*y\\ 2 , oo + 2 I \\C* (y(s) - y h (s)) \\ds. 



As C*C G C(V (A 2 )) DC(V (A )) and \\z\\ 2;O0 
last relation becomes 



tt^-(0)-z7(0)|| <Mh»(T + T 2 )\\z \\ 2 + 2 / \\C*(y(s)-y h (s))\\ds. 



\z \\ 2 (since iA is skew- adjoint), the 
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Substituting the above relation and (2.23) in (2.22), we get 



S 3 <MN h [h e (l + r + r 2 )\\z h + J \\C* (y(s) - y h (s)) \\ds ) . (2.25) 



Substituting (2.20), (2.21) and (2.25) in (2.19), we get for all h e (0,h*) 



\zq - z ,h\\ < M 



(jZT +h 6 [l + (l + T + r 2 )N h + tNI]\ \\zo\\ 2 



+ N h / \\C*(y(s)-y h (s))\\ds 
Jo 



which leads to the result (with possibly reducing the value of h*). 



2.2 Full Discretization 

2.2.1 Statement of the main result 



In order to approximate (2.6), we use a finite difference scheme in time combined with 



the previous Galerkin approximation in space. In others words, we discretize the time 
interval [0, r] using a time step At > 0. We obtain a discretization tk = kAt, where 
< k < K and where we assumed, without loss of generality, that r = KAt. Given a 
continuously differentiate function of time /, we approximate its derivative at time tk 
by the formula 

f (t fc ) ~ D t f(t k ) := — . 

We suppose that go,/i G and Ffi, for < k < K, are given approximations of go and 
F(tk) in the space X. We define (g*), for < k < K, as the solution of the following 
problem: for all iph G X h : 



(D t q k h ,<p h ) = ±i(ql<p h ) k - (C*Cq k h ,^ h ) + (F k ,ip h ) 



Th = Qo,h- 



(2.26) 



Note that the above procedure leads to a natural approximation T^ At k of the continuous 
semigroup by setting 



Tjg ^ T+At iJfc go := q k h , 



where q\ solves (2.26) with F k = for all < k < K and for q o h = Tth%- Obviously, this 
also leads to an approximation of L T = T^T^T by setting 



Assume that for all < k < K, y\ is a given approximation of y(tk) in Y and let 
(zfr) and (z^) be respectively the approximations of (2.2) and (2.3) obtained via (2.26) 
as follows: 
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For all < k < K, (z^) = q\ where q\ solves ( |2.26[ ) with = C*y\ and g° = 



For all < k < K, (%) = g£ _/c where q k h solves ( |2.26| ) with F h fe = C*y^ and 



Then, our main result (which is the fully discrete counterpart of Theorem 2.1) reads as 
follows 

Theorem 2.5. Let Aq : T> (Aq) — > X be a strictly positive self-adjoint operator and 
C G C(X,Y) such thatC*C G C (V (A 2 )) n C (V (A )) . We assume that the pair (iA ,C) 
is exactly observable in time r > 0. Let z G T> (Afy be the initial value of ( 2.1[ ). With 



the above notation, let z 0t h,At be defined by (1.7) and denote r] := ||L. 



j t\\C(X) 



< 1. Then 



there exist M > 0, h* > and At* > such that for all h G (0, h*) and all At G (0, At*) 
we have 



ko - Zo,h,At\\ < M 



\ + (h e + At)(l + r)Nl At )\\zoh 

1 — r\ 



K 



+ N hAt AtJ2\\C*(y(t e )-y e h )\ 



E=0 



Corollary 2.6. Under the assumptions of Theorem \2. 5[ we set 

ln{h e + At) 



N, 



h,At 



In t] 



Then, there exist M T > 0, h* > and At* > such that for all h G (0,h*) and 
At G (0, At*) 



\Z0 ~ Z 0) hAt\\ < M T 



(h e + At)\n\h e + At)\\z \\ 2 



K 



|ln(/i e + At)|At^||C*(y(t,)-^)| 



1=0 



. (2.27) 



Remark 2.2. Contrarily to the semi-discrete case, we have not been able to extend our 
results for zq in a larger space than T> (A^). 



2.2.2 Proof of Theorem [231 



The proof of Theorem 2.5 goes along the same lines as the one of Theorem 2.1 in the semi- 
discrete case and uses energy estimates similar to those developed in Fujita and Suzuki 
[HI p. 865]. The main ingredient for the convergence analysis is the following result (the 



space and time) of system (2.6) by (2.26). 



counterpart of Proposition 2.3) which gives the error estimate for the approximation (in 



Proposition 2.7. Given initial states qo G T> (Aq) and qo^ G Xh, let q and q\, for 

< k < K, be respectively the solutions of ( 2.6[ ) and (2.26). Assume that C*C G 
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C (V (A )). Then, there exist M > 0, h* > and At* > such that for all he (0,h*), 
all At G (0, At*) and all < k < K: 



r k 

\Mt k ) - q k h \\ < \\7T h q - q , h \\ + M I At £ \\F(U) - ^11 

{h e + At) [t fc (||gol| 2 +l|F||i,oo + ||F||oo)+^||F|| 



+ 



2,oo 



Proof. Let ri(ifc) denote the residual term in the first order Taylor expansion of q around 
tfc_i, so that 



?(*fc) = ^ aT™ = Dtq ( tk > ~ A~t^ k ' 



(2.28) 



Subtracting (2.26) from the continuous weak formulation (2.6) applied for t — t k and for 



an arbitrary test function (p = iph G X^, we immediately get by using (2.28) that for all 
1 < k < K 

(D t (q(t k ) - q k h ) , <p h ) = ±i (ir h q(t k ), <p h ) k - (C*C (q(t k ) - q k h ) , <p h ) 



+ ^ t (r l (t k ), ( p h ) + {F(t k )-F>:,ip h ). 



The above relation implies that 



[D t (TC h q(t k ) - q k h ) , (f h ) = (D t (TT h q(t k ) - q{t k )) , <p h ) 

± i (rr h q{t k ) ~ q k h , (p h ) k -(C*C (q(t k ) - q k h ) , (p h ) 



+ Xt (ri(t k ), ip h ) + (F(t k ) - F k , <p h ) . (2.29) 



Now, for all 1 < k < K, let 
Using the identity 

1 ,„ „o ,, 



4 = ^hhq(t k )-q h h \\ 2 - 



\u\\ — \\v\\ + \\u — v\ 



2 ) = Re (u — v, u) , Viz, v G X, 



one easily obtains that for all 1 < k < K 

D t £ k < Re (D t (n h q(t k ) - q k ) ,n h q(t k ) - q k ) . 



Substituting (2.29) with ip h = n h q(t k ) — q\ in the above inequality and using the bound- 
edness of C, we obtain the existence of M > such that for all 1 < k < K 

DtSt < [|| A (n h q(t k ) - q{t k )) || + M\\n h q(t k ) - q(t k )\\ 

+ ^l|ri(t*)ll + IW**) " Jtfll] \Whl(h) ~ q k h l (2.30) 
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Using the straightforward relations 



ck-1 



and 



\n h q(t k ) - < V2 + 



(2.31; 



(2.32) 



we obtain from (2.5) and (2.30) that for all h 6 (0, h 



< ML 6 (||Aff(^)||| + ||<K**) Il§) + ^llnftOII + \\F(t k )-Ft\ 



By (2.28) and relations (4.2) and (4.3) in Lemma 4.2 of the Appendix, the last estimate 
yields 



Dtjeji < M^h 6 (1Mb + hWFh^ + ||F|| liC 



+ \\F(tu) - F k h \\ + ^rlMfOlU + h\n{t k )\\ \- (2-33) 



At 



2 AV 



To conclude, it remains to bound the two last terms in the above estimate. By definition 
of 7*1, we have 

ri(t k ) = q(t k -i) - q{t k ) + Atq(t k ), 
in T> {^Aq \ , and thus by the mean value theorem, we get 

lkx(*fc)lli < At SU P U(s)\\i+^t\\q(t k )\\i. 
se[*fc-i,*h] 



Using once again (4.3), we obtain that there exists M > such that 
||?"i(£fc)IU - MAt (\\q0W2 + t k \\F\\ 2 ,00 || ^ || 1,00) • 



(2.34) 



Now by the regularity of q (see Lemma 4.2), the residual 7*1 can be expressed via the 
integral 



in X, and thus 



n(t k ) = / q(s) (Vi - s) ds, 

||n(i & )||<Ai 2 sup ||$( a )||. 

se[tfe_i,tft] 



Using equation (2.4) verified by q and the boundedness of C, we have 

dq 



W)\\= f t (t) 



< 



j t {±iA Q q{t)-C*Cq{t) + F{t) 
?(*)||i + M||g 



Hence, once again by (4.3), we get 

lki(* fc )|| < At 2 f||g || 2 + h\\F\\ 2t00 + ||F|| lj00 + ||F| 



(2.35) 
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Substituting inequalities (2.34) and (2.35) in relation (2.33) provides estimates for D t y^S^ 



I ck 
°h 



£. 



k-1 



At 

(since \\ir h q(t k 



for k = 1, . . . , K, that can be added together to get the desired inequality 



Using this Proposition, we can derive an error estimate for the semigroup T t (for all 



1 < k < K) and for the operator L T = T r (the counterpart of Proposition 2.4) 



Proposition 2.8. Under the assumptions of Proposition \2. 9f the following assertions 
hold true 

1. There exist M > 0, h* > and At* > such that for all h G (0,h*), all At G 
(0, At*) and all < k < K 



||^T+ 9 o-T+ At)fe go|| < Mt k (h 9 + At)\\q \\ 2 . 
W h T7qo-T7 At , qo \\ < M(T-t k )(h + At)\\q o h. 



(2.36) 
(2.37) 



2. There exist M > 0, h* > and At* > such that for all n G N, all h G (0, h*), all 
At G (0, At*) and all < k < K 



(K - K,M,M\ <M[h e + nr (h e + At)} \\q \\ 2 . 



(2.38) 



Proof. 



1. It suffices to apply Proposition 2.7 with F(t k ) = = for all < k < K and 



2. First, we note that 



(2.39) 



Using (2.5), the fact that ||L™||£(x)(A)) < 1 (proved in Lemma 4.1 of the Appendix), the 
first term in the above relation can be estimated as follows 



|L"g -7r,L"g || < Mh e \\q \\ 2 , V/i G (0,h*). 



(2.40) 



For the second term in (2.39), we prove by induction that for all n G N, all h G (0, h*) 

(2.4i; 



and all At G (0, At*) (for some At* > 0) 
By definition, we have 

HtT/iL^O — L^At,fe?o|| 



< 



KhT tk Tlq - T h At)fc T+ g ||, 
(n h T tk -T hAt>k ) ir h T+q \\ 

+ \\ T h,At,k (*hTt k -Tj At fc ) g 1 



Using (2.37) and Lemma 4.1, we get 
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Obviously \\T h At k \\c(x) is uniformly bounded (with respect to h and At), and thus again 
by (2.36) we have 



I 



h,At,k 



(ttaTJ - T+ tk ) g || < Mt k (h 6 + At) \\q \\ 2 . 



So, by adding the two last inequalities, we obtain that 

\\n h l^ k qo - h h>At>k q \\ < Mr (h e + At) \\q \\ 2 , 



showing that (2.41) holds for n — 1. Suppose now that for some n>2 



KL^o - KAMI < M ( n - l ) T i h<) + At ) \\Qo\W 



(2.42) 
(2.43) 



Writing 



n-1. 



we get by using Lemma 4.1, the uniform boundedness of HL^At.fclUpr) with respect to h 
and At, fl2T42~| and ( |2.43[ ) that 

||tt,L™ g - LJ^goll < M [(1 + (n - l))r (// + At)] ||g || 2 , 
which is exactly ( |2.41[ ). Substituting Q2.40D and ( |2.41| ) in ( |2.39[ ), we obtain the result. ■ 



We are now able to prove Theorem 2.5 



of Theorem 2.5 We first introduce the term At K z (0) to rewrite the approxima- 



tion error Zq — Zo,h,At i n the following form: 



n=0 



N h,. 



ZQ — Zq 



hAt =^L>-(0)- ^LJ^fe)' 



n=0 



n=0 



£ L?z"(0)+ J] (L"-L^)^(0) 



?1=0 



n=0 



Therefore, we have 
where we have set 



|| 2o - z ,h,At\\ < S! + S 2 + S 3 , 

n>N h ^ At 
N h , A t 

s 2 =J2 \\(K-KAt,K)z-(o)\\, 



(2.44) 



ra=0 
f N h A t 



So 



EII L 



h,At,K\ 



n=0 



^(o)-fe) ( 
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Since rj = \\h T \\ C (x) < 1, the first term can be estimated very easily 



Si < M 



N h , At + l 



1 — T] 



\Zq 2- 



(2.45) 



The second term S 2 can be estimated using the estimate (2.38) from Proposition 2.8 



( N h, 



S 2 <m\ (h e + nr(h 6 + At)) [||z~(0)|| 2 , Vh E (0, h*),At E (0, At*). 

I n=0 



Therefore, using (1.4), the fact that ||L r ||x)(A 2 ) < 1 (see Lemma 4.1) in the above relation, 
we get that for all h E (0, h*) and At E (0, At*) 



S 2 <M 



1 + (1 + r)N hAt + (1 + T)N 2 hAt ] (h e + At) \\zq\ 



(2.46) 



It remains to estimate the term 5*3. As for the semi-discrete case, on can easily show that 
||I^ft,At,i<:|U(x) is uniformly bounded by 1 (with respect to h and At), and thus we have 



S3 < MN hAt \\z-(0)-( 



< 



MN hAt (\\z-(0) - 7r h z-(0)\\ + ||vr^-(0) - (^ )°||) . 



By using (2.5) and (1.4), we immediately obtain that 

11^(0) -717,2" (0)11 < Mh e \\z \\ 2 . 



(2.47) 



(2.48) 



To estimate the second term 7ThZ~(0) — {z^) , we apply twice Proposition 2.7 first for 
the time reversed backward observer z~{r — ■) and then for the forward observer z + (the 
time reversal step is introduced simply because Proposition 2.7 is written for initial (and 
not final) value Cauchy problems). After straightforward calculation we obtain that for 
all h E (0,h*) and all At E (0, At*) 

||7T & 2-(0) - (^)°|| < M(h 6 + At) [r(||z+(r)|| 2 + \\C*y\\ 1>00 + \\C*y\U 

K K 

+ r 2 ||C*y|| 2i0O ] + At \\C* {y(r - t e ) - y^) || + At £ \\C* (y{t t ) - y{) \\ . (2.49) 

1=1 1=1 

Applying (4.2) of Lemma 4.2 of the Appendix with zero initial data, we obtain that 

|k + (r)|| 2 <T||C*y|| 2)0O . 



As C*C EC(V (A 2 )) nC(V (Aq)) and \\z\\ 2)OQ 
also reads 



\z \\ 2 (since iA is skew-adjoint), (2.49) 



K 



\n h z~(0) - (z-)°\\ < M(h e + At)(r + r 2 )||2 || 2 + 2 At ^ \\C* (y(t e ) - y{) ||. 



1=0 



Substituting the above relation and (2.48) in (2.47), we get 

K 



S 3 < MN hAt I (h e + At)(l + r + r 2 )\\z \\ 2 + At \\C* (y(U) ~ y{) II \ ■ (2.50) 
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Substituting 02145) , 02146) and 02150) in 02T44) , we get for all ft G (0, /i*) and all At G 
(0, At*) 



A' 



P0 — ^0,/i 



Atll < m j N hj & t At J2 \\C* (y(k) - y[) || + \ 



\ z o\\2 



+{h e + At) [l + (1 + r + r 2 )iV hiAt + (1 + r)iVj At 
which leads to the result (with possibly reducing the value of h* and At*). 



\ z o\\2 



3 The wave equation 

Let if be a Hilbert space endowed with the inner product (•, •). The corresponding norm 
of H is denoted by || • ||. Let A : T> (A ) — > H be a strictly positive self-adjoint operator 
and Co G C(H,Y) a bounded observation operator, where Y is an other Hilbert space. 
The norm in V(Aq) will be denoted by || • || a . Given r > 0, we deal with the general 
wave type system 

' w(t) + A w(t) =0, Vt ^ 0, 
y(t) = C w(t), VtG[0,r], 

and we want to reconstruct the initial value (w , wi) = (w(0),w(0)) of (3.1) knowing y(t) 
for t G [0, t]. In order to use the general iterative algorithm described in the introduction, 



(3.1) 



we first rewrite (3.1) as a first order system of the form (1.1). To achieve this, it suffices 



to introduce the following notation: 



z(t) 



A 



w{t) 
w(t) 

/ 

A 
C G C(X,Y) 



X = V [ A§ ) x H, 



V(A)=V(A )xV[Al , . 
C = [0 C ] 



(3.2) 
(3.3) 



The space X is endowed with the norm 



\? II 2 
Pi i 



G A. 



Note that the operator iA is selfadjoint but has no sign so that the problem studied here 
does not fit into the framework of Section [2j We assume that the pair (A, C) is exactly 
observable in time r > 0. Thus, according to Liu [131 Theorem 2.3.], A + = A — C*C 
(resp. A~ = —A — C*C) is the generator of an exponentially stable Co-semigroup T + 
(resp. T~). We set as usually 

L T = t;t+. 

Throughout this section we always assume that (w ,Wi) G V (A 2 ) = V (Ad) x ^(Ad)- 
Thus by applying Theorem 4.1.6 of Tucsnak and Weiss |18j . we have 

w g c {[o,t],v (4)) n c 1 ([o,t],v(a )) n c 2 {[o,r],v (aI 
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The forward and backward observers (1.2) and (1.3) read then as follows (as second-order 
systems) 

w+(t) + A w + (t) + C*C w+(t) = C*y(t), Vt G [0,r], 
w +(0)=0, u>+(0) = 0, 



w~{t) + A w-(t) - C* C w-(t) = -C* y(t), Vt G [0,r], 
w~(r) = w + (t), w~ (r) = w + (r) . 



(3.4) 
(3.5) 



Clearly, the above two systems can be written as a general initial value Cauchy problem 
of the same form (simply by using a time reversal for the second system) 



p(t) + A oP (t) + C*C p(t) = f(t), Vt G [0, t] 
p(0) = po, p(0) = pi 



(3.6) 



where we have set 



for the forward observer ( pL4| ) : /(£) = CqS/W = CqCqw^) and (po,Pi) = (0,0), 

for the backward observer (3.5) : f(t) = 
(po,Pi) = (^ + (r), -w + (r)) G P (A 2 ) = V (4) x P (A 



-Cgj/(r - t) = -C *C w(r - *) and 



Let us emphasize that with these notation, the semigroups T are given by the relations 



Po 




-p(t) 






_p(t)_ 



Po 




" p(r-t) ' 






_-p(r - 1)_ 



(3.7) 



where p solves (3.6) with / = 0. 



In the next two subsections, we propose a convergence analysis of semi-discretized and 



fully discretized approximation schemes for the forward and backward observers (3.4) and 



(3.5 ). Our proof is based on the convergence analysis of the semi and fully discretizations 



of (3.6). As far as we know, the existing literature on the convergence analysis of full 
discretizations of wave-type systems concern only the particular cases of conservative 
systems (i.e. without damping), see e.g. Raviart and Thomas [HI p. 197] or Dautray 
and Lions (U p. 921] and systems with constant damping coefficients Geveci and Kok 
[7] . For a recent review of numerical approximation issues related to the control and the 
observation of waves, we refer the reader to the review paper of Zuazua [19] . 



3.1 Space Semi-Discretization 
3.1.1 Statement of the main result 



We use a Galerkin method to approximate system (3.6 ). More precisely, consider a family 
(Hh)h>o of finite-dimensional subspaces of T> ^4q^ endowed with the norm in H. We 

denote tc^ the orthogonal projection from T> C^oj onto H^. We assume that there exist 
M > 0, 9 > and h* > such that we have for all h G (0, h*) 



177/^ — Lp\\ < Mh 6 \\lf\\ l 



(3i 
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Given (po,pi) G T> {A 2 ), the variational formulation of (3.6) reads for all t G [0, r] and all 
if G T> (AqJ as follows 



(p(f) >¥ >) + (p(t),^)i + (csc p(t),<p) = {f(t),<p) 
p(9)=po, p(o)=pi- 



Vt G [0, r], 



(3.9) 



Suppose that (po,h,Pi,h) G if/, x if^ and //, are given approximations of (poPi) an d / 
respectively in the spaces X and L 1 ([0,r], if). We define as the solution of the 

variational problem 



(p h {t),(p h ) + (ph(t),<p h )i + (C£C p h (t),cp h ) = (f h (t),<p h ) 
Ph(0)=Po,h, Ph(0)=Pi,h- 



Vt G [0, r] 



(3.10) 



for all t G [0, r] and all <^/i G Hh- 

The above approximation procedure leads in particular to the definition of the semi- 
discretized versions T± of the semigroups that we will use. Indeed, we simply set 



T 



h.t 



Po 

Pi 





Ph(t) 




Po 




Ph(r - t) 




Ph(t)_ 


Pk 




_-Ph(r - t)_ 



(3.11) 

The semi-discretized 



where p h solves ( |3.10| > for f h = and (p ,h,Pi,h) = WhPo^hPx) 
counterpart of L T = T^T^ is then given by 

W = T£ T T+ T . 

Assume that is an approximation of the output y in -^ 1 ([0, r], y) and let u;^" and 



denote the Galerkin approximations of the solutions of systems (3.4) and (3.5), satisfying 
for all t G [0, r] and all <fh G fifc. 



w+(0) = 0, <(0)=0, 2 

[w^(t),ip h ) + (w h {t), l f h ) 1 _ - (C%C Wk(t),<p h ) = - (C*y h (t),p h ) , 
v h( T ) = w h( T ), *%(r) =w£(t). 
With the above notation, the main result of this section reads as follows. 



(3.12) 



(3.13) 



Theorem 3.1. Let Aq : T> (Aq) — >■ H be a strictly positive self-adjoint operator and 



C G £(H,Y) such that C^C G £ [V L4| J J n £ YD [A$ J J. £>e/me (A,C) 5?/ (3.2) 



and ( 3.3[ ). Assume that the pair (A,C) is exactly observable in time r > and set 
7] := ||L r ||£(x) < 1. Let (w ,Wi) G Z> (^o) x ^ (A)) ^ e ^ e initial value of (3.1) and Zei 
(w 0>h ,w lth ) be defined by 



Wo,h 

Wl,h 



n=0 



(3.14) 
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Then there exist M > and h* > such that for all h G (0, h* 



2 



1 — T] 



+ h e rNl ( \\w \\a. + ||wi||i 



+ iV h / ||C *(y(s)-i/ h (a))||ds 
'o 



Corollary 3.2. Under the assumptions of Theorem 3.1, we set 



In 7/ 



Then, there exist M T > and h* > such that for all h G (0, h*) 



\wq - w ,h\\i + \\wi - w lth \\ < M T 



h e In 2 h ^||wo|| s + ||wi 



+ \\nh\ \\C*(y(s)-y h (s))\\ds 
'o 



(3.15) 



3.1.2 Proof of Theorem PI 



The next Proposition provides the error estimate for the approximation of (3.9) by using 



the Galerkin scheme (3.10) 



Proposition 3.3. Given (po,Pi) G V [AqJ x V (A ) and (po,h,Pi,h) H h x H h , let p and 

Ph be the solutions of §3} and ^ecfciy. 04 C 'C„ 6 C (v (a\ 

Then, there exist M > and h* > suc/i t/iat /or a// £ G [0, r] and a// /i G (0, h*) 

lkfcp(*) -Ph(*)IU + lkfcp(*) -Ph(*)ll < M i IkfcPo -poaIU + IkfcPi -pi,/ill 



MlbolU + lbxlli + H/IU^ )+i 



l,oo 



+ / \\f(s) - f h (s)\\ds. 



Proof. First, we substract (3.10) from (3.9) to obtain (we omit the time dependence for 
the sake of clarity) for all ip h G H h 

(p ~ Ph, <Ph) + (P ~ Ph, <ph) i + (Co Co (p - Ph) , <Ph) = (f ~ fh, <Ph) ■ 
Noting that (n h p — p, tp h ) i = for all <ph G H h and that Tihp makes sense by the regularity 

2 

of p (this is a direct consequence of relation (4.1) from Lemma 4.2 used with q 
we obtain from the above equality that for all iph £ Hh 



(K h p - Ph, Ph) + (nhP - Ph, <fh) 1 = (nhP ~ P, <Ph) + (Cq C (p h - P) , Ph) + (f ~ fh, ¥h) ■ 

(3.16) 
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On the other hand, setting 



we have 



!„ „2 !„ „2 

£a = - 7T h p + -\\iThP -Pa I- 

Z Z 2 



£a = (ttaP - Pa, KhP - Pa) + (ttaP - Pa, ttaP - p h ) i 



Applying (3.16) with <ph = ithP — Ph and substituting the result in the above relation, we 



obtain by using Cauchy-Schwarz inequality and the boundedness of Co that there exists 
M > such that 

£h < (IKap-pII + M\\K h p + ||/ -AH) \\-n-hP - Ph\\ ■ 



Since 



d 



dt 



y2£h, the integration of the above inequality from to t yields 



Whp{t) -pa(*)IU + IKapXO -pa(*)II < M { IkftPo -po,aIU + IKapi -pi,aI 



+ / (\\n h p(s) -p(s) || + \\n h p(s) -p( s )\\)ds+ / ||/( a ) - f h (s)\\d s y (3.17) 



Thus, it remains to bound \\nhp(t) — p(t)\\ and \\7ihp(t) — p(t)\\ for all t G [0,r]. Using 



(3.8) and the classical continuous embedding from V(A a ) to V(A^) for a > (3, we get 



that 



\7T h p(t)-p(t)\\ <Mh e \\p(t)\\l, 

\*hP(t) -p(t)\\ < Mh e \\p(t)\\± < Mh e 



Vi G [0, r], lie (0,/i*). 
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P 



and the right-hand side F 



Using relations (|4.3|) proved in Lemma 4.2 of the Appendix for the first order unknown 
P 



we get for all t G [0, r] and all h G (0, h*) 



WhP(t) -P(t)\\ + hhPit) -p(t)\\ < Mh" \\p \\s + llpiHl +t||/||l,oc + 



Substituting the above inequality in (3.17), we get the result. 



Thanks to the last result, we are now in position to derive an error approximation 
for the semigroups T^ 1 and for the operator L t = TjTf . This result has been recently 
proved in the preprint [2] but we prefer to include the proof for the sake of completeness 
and clarity. 



Proposition 3.4. Let Uh = 

following assertions hold true 



7TA 0' 

n h 



. Under the assumptions of Proposition 



3.3, the 
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1. There exist M > and h* > such that for all t E (0, r) and all h E (0, h*) 



(n fc T7 - t- , 



Po 
Pi 



Po 

Pi 



<Mf/i e (Jbo||§ + IN|i 
<M(r (||pb||| + ||Pi||i) ■ 



(3.18) 
(3.19) 



2. There exist M > and h* > suc/i i/iai /or a// n E N, all t E [0, r] and a// 
fa G (0, h*), we have 



Po 

Pi 



<M(l + nr)/i» (Hpbll. + IMi). 



(3.20) 



Proof. 

1. From relations (3.7) and (3.11) defining the continuous and the semi-discretized 
semigroups and T^, (3.18) and (3.19) follow immediately from Proposition 



2. Let t E [0, r] and (p ,Pi) eV(Aq) xV (A ). We have 



3.3 



Po 
Pi 



< 



:k - iw) 



Po 

Pi 



(IW - K t ) 



Po 

Pi 



(3.21) 



Using (3.8) and the fact that ||Lt||£(x>(A)) < 1 proved in Lemma 4.1 of the Appendix, the 
first term in the above relation can be estimated as follows 



Po 

Pi 



<Mh° Ibolli + lbiHi , Vhe(0,h*). 



(3.22) 



For the second term in (|3.21|), let us prove by induction that for all n E N 

W - Kt) 



Po 

Pi 



<MnrK> IbolU + lb 



i i 



V/i G (0,/i*),Vt e [0,r]. (3.23) 



We have 

(n fe L t - l m ) 



Po 

Pi 



< 



(n,T-T+ - T M T+ 

(n /t T- - Tr t )Tj 



Po 









Po 

Pi 



Tr t (T+-T+) 



Po 
Pi 



By Lemma 4.1 of the Appendix and equation (3.19), we get 



(L\T- - T- )T+ 



Po 









< M(r -t)^(||Po||| + lbi||i), Vhe(0,h*). 



Obviously ||T ft \\c(x) is uniformly bounded (with respect to h), and thus by (3.8) and 
equations (3.18)-(3.19), we have (by using the triangle inequality) 



Po 
Pi 



<M^(||p ||| + lbi||i), Vhe(0,h*). 
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Consequently 

(U h h t -h h ,t) P ° < Mrh e (\\po\\ a Id), V/iG (0,/i*),Vt G [0,r] (3.24) 

which shows that (3.23) holds for n = 1. Suppose now that for a given n > 2, there holds 

< M(n - l)rh e (|bo||§ + lbi||i) • (3.25) 



(W 1 - l»7 14 







Pi. 





From 







< 


Pi. 













+ 




Pi. 











Pi. 





we obtain by using (3.24) and (3.25) that (thanks to Lemma 4.1 and the uniform bound- 
edness of HL^tH^x) with respect to h) 



(II fc L? - L? , 



Po 




Pi. 





< Mnr/i e (j|po||§ + |bi||i 



which is exactly (3.23). Substituting (3.22) and (3.23) in (3.21), we obtain the result. 



Now, we can turn to the proof of Theorem 3.1 



N h 



of Theorem 3.1 Introducing the term ^^L^ T 



n=0 



'ur(o) 

w (0) 



w 



5>; 

n=0 



'ty-(O)' 
w-(0) 



n=0 



'«fi(0)' 



, we first rewrite the error term 
in the following form 



w () 

Wi 



Wo,h 



n>N h 



'uT(O)' 
iir(0) 



N h 



n=0 



n=0 



'uT(O)' 
ti;-(0) 



't£T(0) -«£(())' 
ur(0) -«£(()) 



Therefore, we have 



where we have set 





Wo 




wo,/. 






3. 




_Wl,h 





<s 1 + s 2 + s 3 , 



(3.26) 



n>N t 
N h 

s 2 =j: 



w-(0j 
tir(O) 



n=0 



'tw-(O)" 
u»-(0) 



5, 



v n=0 



£(X) 
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'iiT(O)' 



Note that the term Si is the truncation error of the tail of the infinite sum (1.5), the 
term S 2 represents the cumulated error due to the approximation of the semigroups T 



w~(0) 
iir(0) 



of the 



while the term S3 comes from the approximation of the first iterate 

reconstruction algorithm. 

Since rj = \\L T \\c(X) < 1> the first term can be estimated very easily using relation 



(1.4): 



Si < M 



rj 



N h + 1 



1 — rj 



\Wq 3 + \\wi\U 



(3.27) 



The term S2 can be estimated using the estimate (3.20) from Proposition 3.4 
S 2 < M h e (||tw-(0)||s + ||uT(0)||iJ , V/16 (0,h*). 

\n=0 J 



Therefore, using (1.4) and the fact that ||L r ||x)(A 2 ) < 1 (see Lemma 4.1) in the above 
relation, we finally get that 



S 2 < M\l + (l + r)N h + N^r\h e (\\w \\* + ||wi||i) , Mh e (0,h* 



(3.28) 



Finally, let us estimate the term S3. As HL^H^x) is uniformly bounded by 1 with respect 
to h, we have 



S 3 < MN h {\\w-(0)-w^(0)\\ h + \\w-(0)-w^(0)\\) 

\ 2 2 

+ ||«r(0) - tt^-(0)|| + \\ic h w-(0) - wl (0)|| J. 



(3.29) 



By using (3.8) and (1.4), we immediately obtain that 



|u>~(0) - 7T ft uT(0)||i + \\w~(0) - ic h w-(p)\\ < Mh e f|KH§ + |K||i) 



(3.30) 



3.3 



To estimate Wtt^w (0) — w h (0)|| 1 + \\rr h w (0) — w h (0)||, we apply twice Proposition 
first for the time reversed backward observer w~(t — ■) and then for the forward observer 



w + (the time reversal is introduced just because Proposition 3.3 can only be applied to 



initial value Cauchy problems). After straightforward calculation we obtain that for all 
h G (Q,h*) 

K«r(o) - ^(o)|| 1 + K«r(o) - ^(0)|| 



< Mh 9 



r(|| W + (r)||a + ||^ + (r)||i + + r^yW^ 



+ 



\C*(y(s)-y h (s))\\ds. (3.31) 



Applying (4.2) of Lemma 4.2 of the Appendix with zero initial data, we obtain that 



k + (r)||| + ||t& + (r)||i<r||C5y||i, e 
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Therefore (3.31) also reads 

Kur(O) - (0)|| | + \\n h w-(0) - w-(0)|| < M//(r + r 2 )||C*y|| liC 



+ 



\C* (y(s)-y h (s))\\ds. 



As CqCo G £ (v (4d)) n£ (l? (^o)) and IMIf.oo + IHkoo = + ||ix;i||i (since A 

is skew-adjoint), the last relation becomes 

\\n h w-{0)-w h {0)\\i + ||7r h tt>-(0)-«>fc(0)|| < Mh 9 {r + T 2 ) (||w ||| + |K||i) 

+ [ T \\C*(y(s)-y h (s))\\ds. 



Using the above relation and (3.30) in (3.29), we get 

S 3 < MN h (h e (l + r + r 2 ) (|K|| f + |H| X ) + jT ||C* (y(s) - y h {s)) \\dsj . (3.32) 
Substituting $27\ , fl3~28j ) and fl3~32| in fl3~26| , we get for all he(0,h*) 



\w - w ,h\\i + \\wi - w lth \\ < 



M 



-,N h +l 



1—7] 



\wo\\ 3 + kui i 

2 



h 6 + + r + r 2 )N h + TN 2 h }^j ( 

+ N h [ T \\C* (y(s)-y h (s))\\ds 



which leads to the result (with possibly reducing the value of h*). 



3.2 Full Discretization 

3.2.1 Statement of the main result 



In order to approximate (3.9) in space and time, we use a finite difference scheme in time 



combined with the previous Galerkin approximation in space. We discretize the time 
interval [0, r] using a time step At > 0. We obtain a discretization t k = kAt, where 
< k < K and where we assumed, without loss of generality, that r = KAt. Given a 
function of time / of class C 2 , we approximate its first and second derivative at time tk 
by 

/(**) - /(t fc _o 



f'(t k ) ~ D t f(t k 
f"(t k ) ~ D tt f(t k ) :-- 



At 

f(t k ) - 2/fa_i) + /fa-g) 
At 2 



We suppose that (po,?i,At;Pi,/i,At) £ -f^/i x -f^/i and for < k < K, are given approxima- 
tions of (po,Pi) and in the space X and if respectively. We define the approximate 
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solution (p^)o<k<K of (3.9) as the solution of the following problem: p\ G H h such that 
for all (p h G H h 

{Duplifh) + (pifh), + {C* Q C D t p k h ,<p h ) = (fl<ph) , 2 < k < K 



P° h = Po,hAt, P l h = Ph + &tp 1AAt . 



(3.33) 



Note that the above procedure leads to a natural approximation T^ At k of the continuous 



operators Tf by setting 



< 


Po 

pk 


~ T+ 

— R h,At,k 


Po 

Pi. 






Po 


— R h,At,k 


Po 
PK 





k ' 

Ph 

DtP k h 



P K h ~ k 
-Dtpf. 



(3.34) 



where p\ solves (3.33) with = for all < k < K and for (poAAt,Pi,h,At) = 
(iihPo, T^hPij- Obviously, this also leads to a fully discretized approximation of the operator 
L r = T~T+ by setting 

^h,At,K = ^h,At,K^h,At,K- 

Assume that for all < k < K, y\ is a given approximation of y{tk) in Y and let (vj^) k 

and (wj l ) k be respectively the approximations of (3.4) and (3.5) obtained via ( 3.33[ ) as 
follows: 

• For all < k < K, (w£) k = p\ where p\ solves (3.33) with ffc = C^y^ and 
(po,hAt,Pi,hAt) = (0,0), 



For all < k < K, («;") 



(p ,hAt,Pi,hAt) = (KT.-AKT) 



p 1 ^ k where p\ solves (3.33) with j ] 



Clvl~ k and 



Then, our main result (the fully discrete counterpart of Theorem 3.1) reads as follows 

Theorem 3.5. Let A : T> (A ) — >■ H be a strictly positive self-adjoint operator and 
C G C(H,Y) such that C* Q C Q G C (t> (MY) nc(v (M)Y Define (A,C) by Q 
and (3.3). Assume that the pair (A, C) is exactly observable in time r > and set 
r] := ||L r ||£(x) < 1- Let (wq,Wi) G T> \AqJ x D (Ao) be the initial value of (3.1) and let 
(w 0A At,w 1A At) be defined by 



W ,h,At 
WihAt 



N h 
n=0 



At,K 



W 



-\o 



Dt{w~ h 



(3.35) 



where D t (w h ) 1 



[w. 



\0 



At 



27 



Then there exist M > 0, h* > and At* > such that for all h G (0, h*) and 
At G (0, At*) 



FO - W ,h,At U + \\ w l - w l,hAt\ 

2 



< M 



V 



N h At + l 



1 — T) 



+ (h e + At) (l + r)Nl M ) (|K||§ + IK| 



A" 



+ iV fciAt At^||Q(^)-^)| 



f=0 



Corollary 3.6. Under the assumptions of Theorem \3.5\ we set 

\n(h 6 + At) 



N, 



h,At 



In rj 



Then, there exist M T > 0, h* > and At* > such that for all h G (0, h*) and 
At G (0, At*) 



I W - W0,/i,At|U + ~ w l,ft,A*ll ^ M r 



(h e + At) \n 2 (h + At) (||w ||§ + |K| 



K 



+ \in(h e + At)\ At J2P*o {yfa) - vi) 



1=0 



(3.36) 



3.2.2 Proof of Theorem ET5l 

As in the semi-discrete case, the main ingredient for the convergence analysis is the 



following result (the counterpart of Proposition 3.3) which gives the error estimate for 



the full approximation of the general system (3.9) by (3.33). 



Proposition 3.7. Given (p ,pi) G V [A^j x V (A ) and (p , h ,At,Pi,hAt) G H h x H h , 
let p and {p^)k be the solutions of (3.9) and (3.33) respectively. Assume that CqCq G 



C [V [Aojj. Then, there exist M > 0, h* > and At* > such that for all 1 < k < K, 
all h G (0, h*) and all At G (0, At*) 

\WhP{tk) -Ph\\\ + \WhP{tk) - DtPhW < M^\\n h p - p 0thAt \\i + \\it h Pi -Pi,h,At\\ 

+ (h° + At) [t k (|bo||a + Hftll! + ||/|| |i0C + ll/IU) +^||/|| 1 ,oo] 

+ At£||/(t,)-^||}. 



Proof. Denote by ri(tk) the residual term in the first order Taylor expansion of p around 
tk-i- Then 

P(tk) = P(tk) ~! (tfc ' l} - 4-Mh) = D tP (t k ) - (3.37) 



At 



At 
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At 



We have 

\\-K h p{h) - D t p k h \\ < \\n h p(t k )-n h D t p(t k )\\ + \\D t (n h p(t k )-p k )\\ 
< ^ t \\ri{t k )\\ + \\D t ^ h p{t k )-p k h )\\ 

Therefore, the error we need to bound satisfies 



hhP{h)-Pl\\i + hhPih) - D t p k h \\ < 2 x /s k h + ^Wnih 



(3.38) 



where we have set for all 1 < k < K 



£ k h = 1 {|| A (n h p(t k ) -p k )\\ 2 + \\-K h p(t k ) -p^||i} . 

On the other hand, if r2(t k ) denote the residual term first order the Taylor expansion of 
p around t k -i, then 



<i{tk) 



At aT 2 ^ 

p(t k ) - 2p(t fe _i) +p(t fc _ 2 ) 

At 2 



(ri(t k ) - ri(t fc _0) - -^^(h), 



yielding 
where 
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p(t k ) = D tt p{t k ) - 7 A 



(3.39) 



Using (3.37) and (3.39), the variational formulation ( |3.9 ) written for t = t k and for an 
arbitrary test function tp — tph £ takes the form 

(D tt p(t k ) -J k ,<fh) + (^hP(tk),(fh)i + (CoC (D t p(tk) - ^ r i(^));V^) = (f(tk),(fh)- 
Subtracting (3.33) from the above relation implies that for all 2 < k < K and all (fh € Hh- 
(D tt (p(t k )-p k ),<p h ) + (7r h p(t k )-p k h ,tp h )± = - (C*C D t (p(t k ) - p k h ) , <p h ) 

+ (l k ^h) + ^ (C* C ori (t k ),<p h ) + (f(t k ) - f k , Vh ) . 
From the above relation, we get that 

(D tt (ir h p(t k ) - p k h ) , <p h ) + (iT h p(t k ) - p k , ip h ) i = (D tt (n h p(t k ) - p{t k )) , <p h ) 



- (C*C D t (p(t k )-fi),cp h ) + < 7 fc ,^> + ^ (C*C Q r x {t k ),y h ) + (f(t k ) - f k , Vh ) . 

(3.40) 

Using the identity 

1 /ii nO ii nO 



■u 



\u 



— v || 2 ) = Re (u — v , m) , Viz, v E H, 
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one easily obtains that for all 2 < k < K 
D t S k < (D tt (n h p(t k )-p k ),D t (ir h p{t k ) -p k h )) 



'^hP{tk) ~Ph,D t (7T h p(t k ) -P k h ))x ■ 



Taking (p^ = D t (nhp(t k ) — p k ) in (3.40) and substituting in the above inequality and 
using the boundedness of Co, we obtain the existence of M > such that for all 2 < k < K 



D t £ k < M 



I At (TTfcpfo) -p{t k ))\\ + || A (ir h p(t k ) -p(h))\\ + h 

1 



At 



\ri(t k )\\ + wm) -m 



D t (7r h p(t k )-pt) . (3.41) 



Using relations (2.31) and (2.32), we obtain from (3.8) and (3.41) that for all h G (0, h 

D t sjs~ k < mL (\\D uP (t k )\U_ + \\D tP (t k )\\ k ) + h k \\ + ^- t \\n(t k )\\ + ||/(t fc ) - f k 



By (3.37), (3.39) and relations (4.2) and (4.3) in Lemma 4.2 of the Appendix for the first 



order formulation of (3.6), the last estimate yields 



Dty/ejl < M\h e (Jbollf + Ibilli + tfcI|/I|x,oo + 



+ \\f(h)-fh\ 



+ 



n(t k ) - ri(t fc _i)||i + — (||ri(t fc )||i + \\r 2 (t k )\\i) 



At 2 ' 



+ 



^||ri(t fc ) - n(t k ^)\\ + ^(lh(4)|| + ||r 2 (t fc )||)|. (3.42) 



To conclude, it remains to bound the terms including the residuals r x and r 2 in the above 
estimate. By definition of r 2 , we have 

r 2 (*fe) = p(tk) - p(t k -i) - Atp(t k ), 
in T> ( Aq ) , and thus by the mean value theorem, we get 



M^aOH i < At sup 



s)\\i+At\\p(t k )\\i. 



Using once again (4.3), we obtain that there exists M > such that 



\r 2 (t k )\\i < MAt (\\po\\z + lbi||i+4||/||i,oo + 



(3.43) 



Now by the regularity of p (see Lemma 4.2 applied to the first order formulation of (3.6)), 
the residual r 2 can be expressed via the integral 



r 2 (t k ) 



d 3 p 



s) {t k -i - s) ds, 
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in if, and thus 



\r 2 (t k )\\ < At 2 sup 

se[t k _ 1: t k ] 



d 3 p 
ds~ 3 ~ { 



s) 



Using equation (|3.6|) verified by p and the boundedness of Co, we have 

dp 



d 3 p 
dt 3 



(t) 



dt 



(t) 



d 
dt 



-A oP (t)-C* C p(t) + f(t) 



<\\p(t)\\i+M\\p(t)\\ + \\f(t)\\. 



Hence, once again by (4.3), we get 



r 2 (t k )\\ < MAt 2 (UpoIU + |bi||i + tfc 



1,00 



For the term implying n, we note that 



ri(t fe ) 



p(s)(t k -i - s)ds, 



tfc-i 



in T> \ Aq j. Hence, by a similar argument and (4.3), 

\\ri(h)\\ < MUri^lli < MAt 2 (\\p \\s + \\ Pl \\i+t k 



1,00 



+ L/i 



Then, we write in D ( Aq 



n(t k ) - ri(t fe _i) 



(p(s - At) - p(s)) (t fc _ 2 - s) ds. 



tk-2 



Using the above relation, it comes by using once again (4.3) 
IM4) - »"i(*k_i)| 



< MAt 



sup 

S6(tfc_2,*fc-l) 



Finally 



< MAt 2 (|bo||| + |bi||i + 



|ri(t fc )-ri(**-i)|| < At 



'tfc_2 Js—At 

< MAt 3 sup 

se(tfc_3,t fe _i) 



d 3 p 
d^ 



da ds, 



Using (3.44) and (4.3), we get 

\\n(t k ) -ri(fjfe_i)|| < MAt 3 (\\p \\s + |bi||i +t fc 



1 J 1,00 



(3.44) 



(3.45) 



(3.46) 



(3.47) 



(3.48) 



Substituting (3.43), (3.45), (3.46), (3.47) and (3.48) in relation (3.42) provides estimates 



for Av^ 



At 



for k = 1,...,K. By adding all these inequalities, we 



immediately get an upper bound for and thus the desired inequality thanks to 

p8| and KM. ■ 



31 



Using this Proposition, we can derive an error estimate for the semigroup T^r (for all 



< k < K) and for the operator L T = T r T+ (the counterpart of Proposition 3.4). 



Proposition 3.8. Let U h = 

following assertions hold true 



n h 
n h 



Under the assumptions of Proposition 



3.7, the 



1. There exist M > 0, h* > and At* > such that for all h E (0,h*), all At E 
(0, At*) and allO<k<K 



(n*T+ - t+ ) 



Po 




Pi. 





(n h T t - - t- 



At,k> 



Po 









< Mt k {h e + At) (||po||| + |bi||i) . (3.49) 
< M(r - t k )(h e + At) (||po||3 + |bi||i) . (3.50) 



2. There exist M > 0, h* > and At* > such that for all n E N, all h E (0, h*), all 
At E (0, At*) and all < k < K 



Po 




Pi. 





< M [h e + nr(h e + At)] (\\po\\> + \\pi\\i ) . 



(3.51) 



Proof. 



1. From relations (3.11) and (3.34) defining the continuous and the fully discretized 



tt and T% L 
2. First, we note that 



operators and Tj^ At k , (3.49) and (3.50) follow immediately from Proposition 



3.7 



J^t k ^h,At,, 



Po 




< 


Pi. 







(L» - U h Ll 



Po 




+ 


Pi. 







Po 




Pi. 





(3.52) 



Using (3.8), the fact that ||L"||£(X)(A)) < 1 (proved in Lemma 4.1 of the Appendix), the 
first term in the above relation can be estimated as follows 



Po 
Pi 



<M^(|MI| + IM|i), Vhe(0,h*). (3.53) 



and all At E (0, At*) (for some At* > 0) 



For the second term in (3.52), we prove by induction that for all n E N, all h E (0,h*) 

(3.54) 



By definition, we have 



Po 









Po 




Pi. 





< Mnt k (h 9 + At) (||po||| + Ibilli) 



< (n fc T t - - Tr )n fc T+ 



Po 




Pi. 





Po 




Pi. 





+ 



T 



h,At,k 



(n h T+ - t+ 



Po 




.■Pi. 
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Using (3.50) and Lemma 4.1 we get 



(n,T- - T- Atik )u h T+ 



Po 




Pi. 





< M(r - t k )(h e + At) (j|po||| + HpxIIx) . 



Obviously HT^ At k \\c(x) is uniformly bounded (with respect to h and At), and thus again 
by (3.49) we have 



So 



(ri/iiL ife — L^At,, 







-■Pi. 





Po 




Pi. 





< Mt k (h e + At) (j|po||§ + |bi||i 



<Mr{h e + At)(\\MU + \\pi\\i 



showing that (|3.54|) holds for n = 1. Suppose now that for some n > 2 

< M(n - l)r(h e + At) ( 



Writing 



IPO || 3 



+ Pi 1 



(3.55) 



(3.56) 



Po 




< 


P 1 . 







Po 




Pi. 





+ ||L hiA4 , fc (n h L™:Vo - L^ t>fc po)|| , 

we get by using Lemma 4.1, the uniform boundedness of HL^At^lUpr) with respect to h 



and At, (3.55) and (3.56) that 



Po 









< M [(1 + (n — l))r(// + At)] (||p Q ||| + HpilU) 



which is exactly (3.54). Substituting (3.53) and (3.54) in (3.52), we obtain the result. 



We are now able to prove Theorem 3.5 



N h , A t 



of Theorem 3.5 Introducing the term L 
in the following form 



hAt,K 



n=0 



, we can rewrite 



w 




W ,h,At 


_Wl_ 




_Wi,hAt_ 



w 

Wi 



w ,hAt 



n=0 



=o 



'ur(O)' 
w-(0) 



N h , A t 



At,K 



n=0 



n>N hAi 



AK") 1 

N h , At 
n=0 

N h>A t 
+ h hAt,K 



«r(0) 
w (0) 



'uT(O)' 



n=0 



w-(0)-(w h )° 
w~(0)-D t (w- h 



— \i 
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Therefore, we have 



\W - WohAtW i + ||wi - Wl,h,At\\ < Si + S 2 + S 3 



(3.57) 



where we have set 



Si 



So 




'tir(O)' 
u>-(0) 



'ur(O)' 
w~(0) 



^3 — I /2 ll^,A*,^l 



C{X) 



«r(0)-(ti£)° ' 

^-(o)-AK) 1 . 



Note that the term 5*1 is the truncation error of the tail of the infinite sum (1.5), the term 
5*2 represents the cumulated error due to the approximation of the semigroups T ± while 

~w-(py 



the term S3 comes from the approximation of the first iterate 



ItT(O) 



of the algorithm. 



Since r\ = \\L T \\c(x) < 1j using relation (1.4), the first term can be estimated very 
easily 



S\< M 



1—7] 



Fop + Fi 1 • 



(3.58) 



The term S2 can be estimated using the estimate (3.51) from Proposition 3.8 : for all 
h G (0,h*) and all At G (0, At*) 



S 2 <M\ ^(1 + nr) (/i e + At) (|K(0)||| + ||tir(0)||i) , 



n=0 



Therefore, using (1.4) and the fact that ||L t ||x)(a 2 ) < 1 (see Lemma 4.1) in the above 
relation, we finally get that for all h G (0, h*) and all At G (0, At*) 



S 2 < M 



l + (l + r)iV Al At + iV^ iAt r] (h e + At) (||w ||§ + 



(3.59) 



It remains to estimate the term 53. By the uniform boundedness of ||L/ l) At,.fir||£(x) by 1 
with respect to h and At*, we have 



S 3 < MN hAt (|K(0) - (w h )°\\ 1 + \\w-(0) - AKH) 

< MNh,At ( \\ w ~ (0) ~~ ^hW^ (0) II 1 + \U h w-(0) - (to,)! 

v 11 11 2 2 

+ \\w-(0) - w"(0)|| + \U h w-(0) - D t (w h Y\ 



(3.60) 



By using (3.8) and (1.4), we immediately obtain that 

\\w~(0) - 7T A «j"(0)||i + \\w~(0) - ic h w-(p)\\ <M(h 9 + At) flKHa + IKHi) . (3.61) 



3.7 



To estimate \\n h w (0) — (w h )°\\i + \\iT h w (0) — D t (w h ) 1 ||, we apply twice Proposition 
first for the time reversed backward observer w~(r — •) and then for the forward observer 
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w + (the time reversal is introduced because Proposition 3.7, as it is formulated, concerns 
initial value Cauchy problems). After straightforward calculation we obtain that for all 
h E (0, h*) and all At G (0, At*) 



Kur(0) - K)°||i + ||7r fc ti»-(0) - AW 11 



< M (h e + At) r(|| W +(r)|| f + ||«,+ (r)|| a + »II|,J + r 2 \\C* y\\ 1>oo 



K 



K 



+ AtJ2\K {y(r -£)- y«- £ ) II + At £ II Co* (y( U ) - y{) ||. (3.62) 



Applying (4.2) of Lemma 4.2 of the Appendix with zero initial data, we obtain that 



\w+{t)\\ % + \\w+{t)\\ x <t\\CI V \\^. 



Therefore (3.62) also reads 



Itt^-(O) - K)°|| i + KuT(O) - D t {w h f\\ < M (h e + At) (r + r 2 )||C*y|| 



l,oo 



K 



2AtJ2\\Co{y(tt)-yl) II- 



1=0 



As C *C ec(v (j$y)nc(v (A*)) and Hlf,oo + \\w\\ hoo = |K||§ + IKIli (since A 
is skew-adjoint), the last relation becomes 



K«r(o) - K)°||i + ||7r fc «»-(0) - AW 11 



< M (// + At) (r + r 2 ) ( 



K 



3 + 



2AtX)l|CS(l/(^)-«J) II- 



£=0 



Substituting the above relation and (3.61) in (3.60), we get 



S 3 < MN hM 



(// + At) (1 + T + T 2 ) (|K||| + IMIi) 

K 

+AtJ2\\Co{y( t *)-yl) II 



(3.63) 



Substituting gjg , ( gjgp and gjg) in g57| ), we get for all /i e (0, h*) and all At e 
(0, At*) 

\\W - Wo,h,At\\i + \\ W l - w l,hAt\\ < 

M f^— -+(h 6 + At) [l + a+r + r^ + r^j^Klla + IKIl!) 

+ iV,, A4 At^||C*(y(t,)-^) || 
e=o 

which leads to the result (with possibly reducing the value of h* and At*). 
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Appendix 

Let A : V (A) 4la skew-adjoint operator and C E C(X, Y) such that C*C E C(V (A)). 
Assume that A — C*C generates a Co-semigroup of contractions on X. 

Lemma 4.1. The operator A — C*C generates a Co-semigroup of contractions on V (A) 
andV(A 2 ). 

Proof. As C E C(X,Y) is bounded, we clearly have V (A) = V(A- C*C). Moreover, 
C*C E CiV (A)) implies that V (A 2 ) = V ((A - C*C) 2 ). The result follows then from 
[TBI Proposition 2.10.4]. ■ 

Lemma 4.2. Given q E V (A 2 ) and F E C ([0, r], V (A 2 )) n C 1 ([0, r], V (A)), let q 
denote the solution of the initial value problem 

f q(t) =Aq{t)-C*Cq(t) + F(t), t E (0,r), 
I q(0) = qo- 
Then, we have the following statements 

1. Regularity: 

qEC([0,r],V {A 2 )) PI C 1 ([0, r],T> (A)) D C 2 ([0, t],X) , (4.1) 

2. Bound for q: 

||?(*)ll«< IkolU + i^lkoo, /or a = 0,1, 2, (4.2) 
5. Bound for q : there exists M > sitc/i £/iat 

||g(*)||« < M(||g |U+i +t||i r || for a = 0, 1, (4.3) 

w/jere ||F|| a)D0 = sup ||F(*)|| a . 

te[o,r] 

Proo/. 

1. By [18, Theorem 4.1.6], we have g G C ([0, r],X> (A 2 )) flC 1 ([0, r],X> (A)). But since 
C*C EC{V (A)) and F E C ([0, t],V (A 2 )) n C 1 ([0, r], D (A)), we have 

(A - C*C) q{t) EC([0,r],V (A)) H C l ([0, r], X) . 

The last inclusion follows then from the fact that q(t) = (A — C*C) q(t) in T> (A). 

2. By Duhamel's formula, we have 

\\q(t)\\ a = T t q + [ T t - s F{s)ds 
Jo 



< l|T t g || a + / \\T t -sF(s)\\ a ds, 

Jo 

< \\q \\ a + t\\F\ 



I 0,00 ) 



where we have used Lemma 4.1 for the last inequality. 

3. Using the estimate (4.2) obtained for q(t) and the continuity of the embeddings 
V (A 2 ) <->V(A)^-> X, we easily get 

\\q(t)L = \\(A-C*C)q(t) + F(t)\\ a , 

< \\q(t)\\ a+1 + M\\q(t)\\ a + ||F|U,oo, 



< M(||g || a+1 + t||F|U +1>oo ) + 
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